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Abstract. - We explore the effects of counterion condensation on fluid- fluid phase separation 
in charged colloidal suspensions. It is found that formation of double layers around the colloidal 
particles stabilizes suspensions against phase separation. Addition of salt, however, produces an 
instability which, in principle, can lead to a fluid-fluid separation. The instability, however, is 
so weak that it should be impossible to observe a fully equilibrated coexistence experimentally. 



Colloidal suspensions present an outstanding challenge to modern theories of statistical 
mechanics. In spite of the work extending all the way to the beginning of previous century, our 
understanding of these complex systems remains far from complete. Even such basic questions 
as what is the form of the interaction potential between two colloidal particles inside a solution, 
still remains controversial Q. The suspensions most often studied experimentally consist of 
polystyrene sulphate spheres with diameter in the range of 10~ 6 — 10~ 4 cm, and 10 3 — 10 4 
ionisable surface groups. The typical solvent is water at room temperaure. The main stumbling 
blocks are the presence of long range Coulomb interactions and the tremendous asymmetry 
existing between the polyions and their counterions. The ratio of the bare charge of a colloidal 
particle to that of its counterion can be as high as 10000:1. This large asymmetry completely 
invalidates most of the methods of liquid state theory, which have proven so successful in the 
studies of simple molecular fluids. 

A particularly interesting question that has provoked much controversy over the last two 
decades concerns itself with the possibility of fluid-fluid phase separation in charged colloidal 
suspensions. A naive argument, based on the theory of simple molecular fluids, suggests that 
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fluid-fluid (or liquid-gas) coexistence is only possible in the presence of sufficiently long ranged 
attractive interactions. Thus, it has been proposed by some authors that a phase separation, or 
even existence of voids in charged colloidal lattices requires attraction between the polyions . 
Although appealing intuitively, this point of view is difficult to justify within the framework of 
statistical mechanics. The fundamental observation is that colloidal suspension is a complex 
fluid for which many body effects play the fundamental role. It is, therefore, erroneous 
to confine attention to pair interactions between the colloidal particles while ignoring the 
significantly larger contributions to the free energy arising from the presence of counterions || . 
This point has also been emphasized by van Roij and Hansen (RH) Q who demonstrated, in 
the context of the linearized density functional theory, existence of "volume" terms, which 
can drive phase separation even for pairwise repulsive interactions S. The prediction of a 
liquid-gas phase separation in an aqueous solution of like-charged colloidal particles seems, 
however, to be contradicted by the recent simulations of Linse and Lobaskin ||, who did 
not find any indication of phase transition in suspensions with monovalent counterions. This 
apparent discrepancy between the simulations and the density functional theory suggests that 
a closer look at the mechanism of phase separation is worth while. Since the first simulations 
were performed in the absence of salt, as a starting point, we shall concentrate our attention 
on this regime. 

Our model consists of N p = p p V spherical polyions of radius a, inside a homogeneous 
medium of volume V and dielectric constant D. Each polyion carries Z ionized groups of 
charge —q uniformly distributed over its surface. A total of ZN p monovalent counterions of 
charge -\-q are present in order to preserve the overall charge neutrality of solution. In the 
absence of salt, the counterions can be treated as point-like. 

All the thermodynamic properties of colloidal suspensions can be determined given the free 
energy. Unfortunately due to complexity of these systems, no exact calculation is possible 
and approximations must be used. We construct the total free energy as a sum of the most 
relevant contributions: electrostatic, entropic, and hard core, / = F/V = f e i + fent + fhc- The 
electrostatic free energy, f e i, is the result of polyion-counterion, f pc , and the polyion-polyion, 
f pp , interactions. Interacions between the monovalent microions are insignificant for aqueous 
solutions and can be ignored Q . 

The polyion-counterion contribution to the total free energy can be obtained in the frame- 
work of Debye-Hiickel theory ||, [?j . Fixing one colloidal particle at the origin, it is possible to 
show that the electrostatic potential in its vicinity satisfies the Hclmholtz equation V 2 tp = K 2 ip, 
where na — (4nZ p*/T*Y^ 2 , and the reduced temperature and density are T* = kgTq 2 D j 'a 
and p* — p p a 3 , respectively. The electrostatic free energy can be obtained from the solution 
of the Helmholtz equation followed by the Debye charging process, yielding 

Z 2 

Pfpc = 2T*(1 + ko) Pp • (1) 

We note that this expression would be identical to the one obtained by RH using the density 
functional theory j|, but for the self energy contribution, flfseif = Z 2 p p /2T*, which wc 
include and RH excluded. Since f se if is proportional to the density, it is irrelevant for any 
thermodynamic calculations, as long as the effects of charge renormalization are neglected. 
Thus, contrary to some earlier claims ||, the Debye-Hiickel theory is fully consistent with the 
density functional theory. 

The polyion-polyion contribution to the free energy is calculated within the variational 
approach proposed by Mansoori and Canfield || . The electrostatic DLVO potential [|l0| is used 
to describe the effective pair interactions between the colloidal particles inside the suspension. 
Based on the Gibbs-Bogoliubov inequality F < F + (U)o, Mansoori and Canfield replace the 
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free energy by the lowest variational bound. The subscript denotes the reference system of 
hard spheres whose diameter plays the role of a variational parameter. The polyion-polyion 
free energy is given by |ll[] 

MM = ^ ! fr|r + ^[A 2 G(A) - 1] . (2) 

The first term is the free energy of the reference hard sphere system, while the second is the 
electrostatic contribution evaluated using the Percus-Yevick pair correlation function for hard 
spheres. The volume fraction, rj, plays the role of a variational parameter. The functions 

G(X) - -VL(A) 



12rj[L(X) + S{X)e x ] ' 
L(A) = i27 7 [(l + |)A+(l + 2r ? )], 

S(X) = (1 - Ty) 2 A 3 + 67/(1 - r?)A 2 + 187? 2 A - 12r?(l + 2rj) , (3) 

are given in terms of A = 2nrQr] 1 ^ 3 , where r$ is a measure of the typical distance between the 
macroions and is given by ^ppr^/i = 1. The variational minimum fj is found by solving the 
equation df pp /dri — 0, this has to be done numerically. The polyion-polyion contribution to 
the total free energy is f PP (fj)- Again we observe that this calculation is very similar to the one 
performed by RH using the density functional theory. We also note that the Eq. (g) already 
includes the background subtractions emphasized by Warren Q . 

The contribution to the total free energy arising from the hard core repulsion between the 
colloidal particles can be approximated by the Carnahan-Starling form p3[, 



P (4 - 30 p ) 
(1 - 4> P f 



phc=p / p ! ) , (4) 



where <j) p = 47r,o*/3 is the volume fraction of polyions. The final contribution to the total free 
energy is due to entropic motion of counterions and polyions and can be expressed using the 
Flory theory Q , 



Pfent = P P In i-jr-j - p p + Zp p In ^ J - Zp p , (5) 

where £ p and £ c are the internal partition functions of the polyions and the counterions, and <f) c 
is the "effective" volume fraction of counterions. Since both the polyions and the counterions 
are rigid, without any internal structure, £ p = 1 and Q c = 1. Although we have assumed 
the counterions to be point-like, the strong electrostatic repulsion will prevent them from 
approaching too close to one another. The distance of characteristic approach d can be obtained 
by comparing the electrostatic and the thermal energies. More specifically from the theory of 
one component plasma we find that d = [(1 + ina/T*) 1 / 3 — 1]/k. The volume fraction 
occupied by the counterions is then (f> c — AirZ p p d 3 /3. Combining all these contributions, we 
obtain the total free energy / = f ent + f pc + f pp + fi lc of a colloidal suspension. We observe 
that when ZXb/cl is bigger than approximately 15, the Hclmholtz free energy is no longer a 
convex function of colloidal density and the separation into two coexisting phases becomes 
thermodynamically favorable. Here A^ = a/T* is the Bjerrum length and is equal to 7.2A 
in water at room temperature. In Fig. [l] we demonstrate some characteristic pressure-density 
isotherms, exhibiting the familiar van der Waals loop, indicating the presence of a first order 
phase transition. 

Is the phase transition found above consistent with the underlying approximations of the 
Debye-Hiickel and the linearized density functional theory? Clearly the fact that the transition 



4 



EUROPHYSICS LETTERS 



5 



2.5 




-2.5 1 ' ' ' ' ' ' ' ' ■ 1 

0.02 0.04 0.06 0.08 0.1 

4» 

Fig. 1. - The pressure-volume fraction, <j> — | 7ra3 ppi isotherms for salt-free colloidal suspensions at 
room temperature, Xb ~ 7.2A. The circles are for Z = 2000, squares for Z — 2200 and triangles for 
Z = 2500. The radius of poly ions is a = 1000 A. Note presence of the van der Waals loop for Z = 2200 
and Z = 2500. 

is located at Z/T* w 15, i.e., the strong coupling regime, should leave us concerned. Certainly, 
in this regime the charge renormalization due to strongly associated counterfoils should play a 
significant role Jl6| . Fortunately, it is fairly straight forward to include the effects of counterion 
condensation directly into the theory presented above ||, [l?]]. To achieve this, we separate 
counterions into condensed and free. For simplicity we shall assume that each polyion has 
an equal number n of condensed (associated) counterions . The density of free counterions 
is then pf = (Z — n)p p . We shall suppose that the only effect of condensed counterions is 
to renormalize the bare charge of colloidal particles, while the screening is performed by the 
free microions. The effective charge of a polyion-counterion complex is then Z e ff = Z — n. 
The free energy f(ri), taking into account the charge renormalization, is obtained by replacing 
Z — > Z — n in all formulas, including the inverse screening length n. In addition, since the 
complexes now have structure ||, their internal partition function can be approximated by 
C P = [Zl/((Z — n)lnl)]e~P E ™ , where E n is the electrostatic energy of n counterions condensed 
onto the surface of a polyion. The electrostatic energy of association can be obtained through 
the charging process j|, yielding f3E n = —{Zn — n 2 /2)/T*. 

For fixed volume and number of particles, the equilibrium state of a colloidal suspension is 
determined by the minimum of Hclmholtz free energy, f ren = min„/(n). Predictions of this 
theory for Z e ff have been recently tested and found to be in a semi-quantitative agreement 
with the experiments p8| . We find that the charge renormalization, Fig. 2b, has a most 
profound effect on the free energy. The f ren is a convex function of colloidal density for all 
values of Z, see Fig. 2a. We conclude that counterion condensation stabilizes salt-free colloidal 
suspension against phase separation. 

We now turn our attention to suspensions in the presence of salt. Since a system containing 
point-like positive and negative particles is intrinsically unstable, we assign to each counterion 
and coion a characteristic radius a c . The concentration of monovalent salt is designated by 
p s . The calculation now proceeds as the one outlined above. The density of free microions is 
pf = (Z — n)p p + 2p s and na = {Airpr/T*) 1 / 2 . The effective charge of a polyion-counterion 
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Fig. 2. - (a) The pressure-volume fraction isotherms for salt-free colloidal suspensions. All parameters 
are as in Fig. lj. The solid curves are for Z = 2000 (circles), 2200 (squares), 2500 (triangles) with 
counterion condensation taken into account. Note that Z e ff is almost the same for all three bare 
charges, and the isotherms collapse onto one monotonically increasing curve. There is no phase 
separation. The dashed curve, included for comparison, neglects the counterion condensation, and 
demonstrates the presence of the van der Waals loop and the phase separation for Z — 2500. (b) The 
effective charge Z e ff vs. Z for various colloidal volume fractions. Note that the effective charge is 
quite insensitive to colloidal concentration. 



complex is Z e ff = Z — n. The polyion-counterion and the polyion-polyion contributions to 
the total free energy are the same as for the case without salt, but with the new definition of 
K, since all free microions contribute to screening. The cntropic free energy is given by 

Pfent = Pp In - Pp+ P+ 1» (jf^ - P++ P- ln (j^j ~ P - ' ( 6 ) 

where the density of free counterions is p+ = (Z — n)p p + p s and the density of coions is 
P- = p s . The volume fractions are (/>+ = 47rp + a^/3 and <f>- = Airp-a^/S- 

To study a possibility of phase separation in this multicomponent system is significantly 
more difficult than for the case of salt-free suspensions. The clearest indication of phase 
transition can still be obtained from the pressure-density isotherms. A caution, however, must 
be taken since the two coexisting phases do not necessarily have the same concentration of salt, 
but must have the same chemical potential. This can be controlled by putting the suspension 
in contact with a hypothetical reservoir containing an aqueous solution of salt. In Fig. ^, 
we present the pressure-density isotherms for suspensions with salt. We see that even an 
extremely low concentration of salt is sufficient to shift the phase transition to the region of 
parameter space where the charge renormalization plays only a marginal role. The suspension 
phase separates, however, the transition is much weaker than it would be in the absence of 
counterion condensation, Fig. |^. The instability region forms a closed loop in the (p s , p p ) plane. 
For the values of Z, a, a c , and As used in Fig. |^ the upper critical point is at ~ 6ApM and 
4>p ~ 0.02. For salt concentrations p s > the suspension is completely stable. The lower 
critical point is located at p l s w 0.28/iM and <fr l p ks 0.000025. Since in practice it is impossible 
to deionise water bellow p s « lpM, all aqueous suspension should — in principle — exhibit 
a miscibility gap at sufficiently low colloidal volume fractions and salt concentrations. From 
Fig. |^, however, it is clear that fluid-fluid transition is so weak that it is highly unlikely that 
a fully equilibrated coexistence can be observed experimentally. 

We conclude that the non- linear effects associated with the presence of double layers strongly 
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Fig. 3. - The pressure-volume fraction isotherms for colloidal suspensions with two different reservoir 
salt concentrations, circles for p a = 25 /iM and squares for p s = 1 /iM; \b = 7.2 A, a — 1000 A, 
a c — 2 A, Z — 2500. The solid curves are when the charge renormalization is neglected, while the 
dashed curves are when the charge renormalization is taken into account. We note that in the presence 
of salt, the counterion condensation does not fully destroy the phase transition, but only makes it 
weaker. Inset amplifies the region of small densities, showing a part of the van der Waals loop for 
p s = 1 p,M. Note the scale of pressure. 



modify the critical phenomena of charged colloidal suspensions. In particular, we find that 
aqueous suspensions with monovalent counterions do not phase separate in the absence of 
salt. Addition of a small amount of 1:1 electrolyte produces a very weak instability. In 
fact the phase transition is so weak that it is unlikely to be observed experimentally. The 
metastable effects associated with the mathematical singularity, however, might indeed have 
been detected, appearing as dilute voids in homogeneously deionised suspensions Q]. 

We would like to stress that our conclusions are only applicable to aqueous suspensions 
with monovalent counterions. In the case of multivalent counterions, the correlations between 
the condensed counterions can become so strong as to lead to an effective attraction between 
the polyions JlSf] , which in turn can drive a liquid-gas phase separation. Similarly, in organic 
solvents of low dielectric constant, the correlations between the condensed counterions can 
become sufficiently strong as to produce phase separation [ pp| . Unfortunately, this correlation 
induced instability cannot be studied at the level of the theory presented in this letter since 
the interactions between the polyions are described by the DLVO potential, which does not 
include correlations. The DLVO potential is sufficient as long as the attention is restricted 
to water with monovalent counterions, but most certainly fails if the electrostatic interactions 
between the condensed counterions are significant, as is the case for low dielectric solvents and 
multivalent counterions. 

The theory presented above suggests that the counterion condensation stabilizes charged 
colloidal suspensions against fluid-fluid phase separation. It, however, says nothing about 
the possibility that suspension freezes, forming an ordered lattice. Fluid-solid transition 
will certainly occur at sufficiently large colloidal volume fractions and will depend on the 
concentration of salt in the system. To study this transition requires an accurate free energy 
of the crystalline state. This can be obtained using the density functional theory within a 
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suitable cell geometry The equality of pressure and chemical potential in the two phases will 
lead to a fluid-solid coexistence curve. 
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